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Abstract 

We discuss how to define and obtain the running coupling of a gauge theory in the 
approach of the Schwinger-Dyson equation, in order to perform a non-perturbative 
study of the theory. For this purpose, we introduce the nonlocally generalized gauge 
fixing into the SD equation, which is used to define the running coupling constant 
(this method is applicable only to a gauge theory). Some advantages and validity of 
this approach are exemplified in QED3. This confirms the slowing down of the rate of 
decrease of the running coupling and the existence of non-trivial infra-red fixed point 
(in the normal phase) of QED3, claimed recently by Aitchison and Mavromatos, 
without so many of their approximations. We also argue that the conventional ap- 
proach is recovered by applying the (inverse) Landau-Khalatnikov transformation to 
the nonlocal gauge result. 
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1 Introduction 



In quantum field theory, a change in the renormalization scale fi accompanied by 
suitable change in coupling and mass does not change the theory. The invariance 
of the theory under such transformation is called the renormalization group (RG) 
invariance. For a given physical theory, we have, for each value of fi, a definite value 
of the coupling and m(/i). These are called the effective (or running) coupling 
and mass. Physical quantities like the S-matrix are invariant under the change of 
variable: (/i, rn(fi)) — > (//, g(fjf), m(//)). This invariance leads to a differential 
equation for the running coupling and mass, called the RG equation Jl||. This is the 
easiest way in practice to compute the effective coupling and mass. The coefficients 
P,j m , ... in the RG equation are called the RG coefficients which express important 
properties of a theory. 

The most important application of the RG is to compute large-momentum, i.e. 
ultraviolet (UV) behavior. However, RG methods can be used to compute infrared 
(IR) behavior too. Certainly this is true in a purely massless theory (or if masses can 
be neglected), and the IR behavior is computable perturbatively if and only if the 
theory is not asymptotic free. But in a massive theory it is not useful to take \i much 
less than a typical mass, for one obtains logarithms of m//i and these prevent simple 
use of the perturbative method when /i <C m, see e.g. [Kl- 
in four spacetime dimensions, it is well known that the only theories that are 
asymptotically free are non-Abelian gauge theories with a small enough number of 
matter fields, see ||. In an asymptotically free theory, the effective coupling g(fi) 
goes to zero when [i goes to infinity, so that short-distance (or high-energy) behavior 
is computable perturbatively. But, when /i is small, g(n) is large, so that IR behavior 
cannot be computed reliably by perturbation theory. This is the case for QCD4. 

If we consider a non-asymptotically free theory in four dimensions, like {4> 4 )a and 
QED 4 , then the effective coupling increases with energy. Thus, in such theories it 
is impossible to compute the true high-energy behavior by perturbation theory. But 
when \i goes to zero, so does the effective coupling. Hence, we can compute IR 
behavior in such a theory, just as we compute the UV behavior in an asymptotically 
free theory. However, note that the coupling in QED is a := e 2 /47r ~ 1/137. This is 
so small that the non-perturbative region in QED4 does not occur until very many 
orders of magnitude beyond experimentally accessible energies. Nevertheless, the 
existence of a non-trivial UV fixed point in a non-asymptotically free theory in four 
dimensions has been suggested by Miransky for strongly coupled QED4 where the 
fixed point is expected to be located at a = a* ~ 0(1) of order unity. In the strong 
coupling region, we need to deal with the theory non-perturbatively. 

As well as simulations based on lattice gauge theory [[|, the Schwinger-Dyson 
(SD) equation has played an important role in the non-perturbative analytical study 
of strongly coupled gauge theories. In the actual analysis of the SD equation, an ap- 
proximation of constant coupling (standing coupling) has been taken together with a 
bare vertex approximation as the simplest approximation, which is usually called the 
quenched planar (or ladder) approximation. However, the introduction of running 
coupling is indispensable to study the unquenched QED, QCD and extended techni- 
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color theory, etc, see e.g. 0. So far, this has been done in most cases in a somewhat 
unsatisfactory way, in my opinion, in which the expression for the running coupling 
constant obtained by (RG-improved) perturbation theory (at most leading logarithm) 
was substituted into the SD equation. However, the running coupling itself should 
be calculated within the SD equation approach. Such a kind of calculation was tried, 
for example, in the massive gauge boson theory in four-dimensions and iV-flavor 
QED 3 m H 0) 0) HH by solving the SD equation for the wavefunction renormaliza- 
tion function A(p 2 ) of the fermion, using either a bare vertex or a reasonably simple 
ansatz for the vertex function. In the original analysis 14] of multi-flavor QED3, 



no wavefunction renormalization was assumed from the beginning, based on a naive 
1/N argument. Quite recently, self-consistent solutions of QED 3 have been studied 
more extensively by solving the coupled SD equation under various ansatzes for the 
vertex Jl5| . 



The choice of vertex ansatz is the most difficult problem in truncating the in- 
finite hierarchy in the SD equation approach. Even if the vertex ansatz might be 
reasonable, we need to make a number of approximations in solving the SD (integral) 
equation, at least analytically. The nature of those approximations is totally different 
from what one encounters in solving a differential equation. Usually, obtaining an an- 
alytical solution for an integral equation is much more difficult than for a differential 
equation. Therefore, the SD integral equation is often converted into a differential 
equation by simplifying the kernel, although the integral equation can not in general 
be transformed into a differential equation in a mathematically rigorous sense. [] 

In this paper we discuss another approach to calculating the running or effective 
coupling in a gauge theory. In this approach we introduce the nonlocal gauge-fixing 
[]16|) |T7| , |18| . Here the existence of gauge degrees of freedom is an essential ingredient, 
so this approach is applicable to a gauge theory only. 

The conventional approach relies on a specific choice of vertex function. Therefore, 
one can not be free from the criticism whether the result is due to an artifact of the 
specific ansatz adopted for the vertex function or not. Of course, it is impossible to 
get rid of this criticism completely. However, we can allow a certain class of ansatz for 
the vertex function which gives a weaker restriction than the conventional approach. 
In this paper, we choose the ansatz for the vertex: 

T,(p,q)= 1 ,G(p 2 ,q 2 ,k 2 ) (1.1) 

where the function G is an arbitrary function of the fermion momenta p 2 , q 2 = (p — k) 2 
and the gauge-boson momentum k 2 . We shall see that the resulting running coupling 
is essentially independent of the explicit form of G, provided the vertex function 
satisfies the Ward-Tahakashi identity. This is a major advantage of this approach. 

So far, the nonlocal gauge has been applied to gauge theories in four dimensions 
17 , TH] , the gauged Thirring model |T^] and QED3 [ 2"0 , |2"T| , |2"2 |. In those works, a 



1 The numerical approach enables us to solve the integral equation as well as the differential 
equation, although the accuracy and the convergence of the algorithm for solving the integral equa- 
tion has not been established so well. Justification of the approximations made in the analytical 
study has been done through the numerical calculations. 
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bare vertex approximation has been taken from the beginning. This assumption is 
not necessary or is weakened as shown in this paper. 

The second advantage of this approach is as follows. This approach is free from 
various approximations which are required to solve the SD equation analytically for 
the wavefunction renormalization, because we never solve it to obtain the running 
coupling constant. Instead, we derive a differential equation which should be sat- 
isfied by the nonlocal gauge function. The nonlocal gauge is obtained by a simple 
quadrature without solving any self-consistent equation. This is an advantage of this 
approach. 

We apply this method to study the non-perturbative behavior of QED3, in partic- 
ular, the IR behavior and the rate of running in the intermediate momentum region 
which have been extensively studied recently by Aitchison and Mavromatos [pi and 



their collaborators [p4}| . We also present the RG-like argument for the behavior of 



the running coupling within this approach. 

This paper is organized as follows. In section 2, we introduce the nonlocal gauge 
fixing and set up the SD equation under a class of ansatzes for the vertex function in 
D(> 2)-dimensional gauge theory. 

In section 3, we derive the differential equation which is obeyed by the nonlocal 
gauge such that there is no wavefunction renormalization A(p 2 ) = 1. Here we in- 
troduce the nonlocal gauge with and without IR (not UV) cutoff. The IR cutoff is 
necessary to discuss the IR behavior of QED3 in the following sections. 

In section 4, we discuss how to obtain the non-trivial wavefunction renormalization 
in the usual (Landau) gauge from the result in the nonlocal gauge. For this, we apply 
the inverse LK transformation. 

In section 5, we define the running coupling through the nonlocal gauge in this 
approach and compare it with that in the conventional approach. Here we see several 
advantages of this approach for obtaining the running coupling. 

In section 6, we study in detail the behavior of the running coupling of QED3 
obtained in the previous section, paying particular attention to the IR behavior. The 
RG-like interpretation of this result is also given here. 

The final section is devoted to conclusion and discussion. 



2 Schwinger-Dyson equation 
2.1 Introduction of the nonlocal gauge 

First, we discuss the SD equation for the fermion propagator. In accord with the bare 
fermion propagator 

S (p) = (Yp,~m )-\ (2.1) 
we write the full fermion propagator as 

S(p) = [A(p 2 )Y% - B(p 2 )]-\ (2.2) 
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For a class of gauge theories in D = d + 1 dimensional space-time, the SD equation 
for the full fermion propagator in momentum space is given by 



S-\p) = S Q \p) + J -^- llx S(q)T u (p, q )D lxl/ (p- q ), (2.3) 

where T u (p, q) is the full vertex function and D^ u (p — q) is the full gauge boson 
propagator. We always use p, q for the fermion momentum and k = p — q for the 
gauge-boson momentum. This class includes QED and the gauged Thirring model 



191, pq] . Note that this SD equation can be decomposed into a pair of integral 
equations for the wavefunction renormalization function A(p 2 ) and the mass function 
B(p 2 ), as shown in the next section. The SD equation for the full fermion propagator 
S(p) should constitute a closed set of equations together with the SD equations for 
the full vertex function T u (p, q) and the full gauge boson propagator D^(p — q) which 
will be specified below. 

Next, in order to specify the gauge boson propagator, we discuss the gauge- fixing. 
In this paper we consider a more general gauge fixing than the usual one, the so- 
called the nonlocal gauge-fixing. In configuration space, the gauge fixing term in the 
nonlocal gauge |2^] is given by 

£gf = ~\f[A{x)\ J d°y_L-jF[^(y)], (2.4) 

with a gauge-fixing function F[A\. In this paper we take the Lorentz-covariant linear 
gauge: 

F[A] = d^Ap. (2.5) 

In momentum representation, the gauge-fixing parameter £ becomes momentum- 
dependent, namely, £ becomes a function of the momentum: £ = £(&). Here it 
should be noted that £ _1 (fc) is the Fourier transform of ^ 1 (x): 

= J ^e^r^), r\k) = J ^e-^r 1 ^), (2.6) 

while £(fc) is not the Fourier transform of £(a;). If £(/c) does not have momentum- 
dependence, i.e., £(&) — > £, then £~ 1 (a; — y) — > S(x — y)^ 1 and hence the nonlocal 
gauge-fixing term (|2.4j ) reduces to the usual gauge-fixing term: 

Cgf = -^(F{A(x)}) 2 . (2.7) 

It is easy to show that the SD equation for the full gauge-boson propagator is 
given by 



D-J(k) = D^-\k)-U, u (k), 

n^(*) := <?f-^trh^(p)r v (p,p-k)S(p-k)], {21 
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where the bare gauge-boson propagator D^J(k) in the nonlocal gauge ( |2.5| ) is given 
by 

D$-\k) = k 2 9fMU - k^K + t(k)-%k v . (2.9) 
In gauge theory, the vacuum polarization tensor should have the transverse form: 

k,,k,. 



n /a/ (*) = Uv-jf Jn(*) 



(2.10) 



provided the gauge invariance is preserved. Hence the full gauge-boson propagator is 
of the form 



D^(k) = D T {k)[ gflv 



ku,k v \ £,(k) k^k v 



k 2 



k 2 k 2 



D T {k) 



1 



k 2 - n(fc) 

In this paper we take another form for the full gauge boson propagator: 



(2.11) 



D^{k)=D T {k 2 ) 



9iw ~ v(k 



k k 



k/j, '■— Pn Q/jt- 



(2.12) 



By comparing (|2.11| ) with (|2.12|) , the correspondence between £ and rj is given in 
momentum space as follows. 



ak) = [i - v(kW - n(£0AT\ v(k 2 ) = i - - n(k)/n (2.13) 



2.2 Vertex ansatz and gauge choice 

Finally, we must specify the full vertex function T u (p, q). Of course, the vertex func- 
tion obeys its own SD equation and all the SD equations should be solved simulta- 
neously. In principle, it might be possible to do that. However, it is rather difficult 
or impossible to actually carry out this scenario, except for some exactly solvable 
models (see e.g.|25|]), due to the infinite hierarchy of the SD equations. Therefore 



we need to truncate the SD equations so that they become tractable analytically or 
numerically. One usually adopt an ansatz for the full vertex function T M , instead of 
solving the SD equation for the vertex function. Such an ansatz is suggested from 
various consistency requirements, see e.g. pUfl . In this paper we adopt the following 



simple ansatz: 

T,( P ,q)=G(p 2 ,q 2 ,k 2 )^, (2.14) 

where g M denotes the momenta of the fermion and k^ the momentum of the gauge 
boson. Here G = G(p 2 ,q 2 ,k 2 ) is an arbitrary function of p 2 ,q 2 and k 2 except for a 
restriction specified below. 

The truncation of the set of SD equations in the gauge theory should be done self- 
consistently in such a way that the resulting truncated SD equations respect gauge 
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invariance (as much as possible). For example, the vertex function should satisfy the 
Ward-Takahashi (WT) identity: 



(p " lY^iP, q) = Sip)' 1 - S(q)-\ (2.15) 

which is a consequence of the gauge invariance of the theory. When the function 
B(p 2 ) vanishes (in the symmetric phase or normal phase) or is extremely small (in 
the neighborhood of the critical point in the broken phase or superconducting phase), 
the WT identity reduces to 

(P " q) ~ [A(pV - A(q 2 )q»} lfl . (2.16) 

This tells us that G(p 2 , q 2 , k 2 ) should be expressed in terms of A(p 2 ), A(q 2 ) and A(k 2 ): 

G(p 2 , q 2 , k 2 ) = F[p 2 , q 2 , k 2 , A(p 2 ), A(q 2 ), A(q 2 )], (2.17) 

Some examples will be shown below. A restriction for our approach to work is 

dG(p 2 ,q 2 ,h 2 ) 



dk 2 



0, (2.18) 



where the derivative is taken with p 2 and q 2 being fixed. This implies that the 
argument of G does not depend explicitly on the angle d of the inner product p ■ q = 
pq cos •& which comes from k 2 := {p — q) 2 = p 2 + q 2 — 2p-q. However, G can depend on 
k 2 implicitly, through A(k 2 ). The reason why we need this restriction will be made 
clear in the next section. In what follows we restrict our discussion to this case. 

The simplest choice is the bare vertex approximation: T M (p, q) = 7 M . This approx- 
imation is usually called the ladder approximation. And it is said that the ladder 
approximation breaks gauge invariance (this is obvious, since Q2.16p can not be sat- 



isfied unless A{p 2 ) = 1 which is impossible as a solution of the SD equation except 
Landau gauge.). This statement is sometimes used as indicating that the ladder ap- 
proximation in the SD equation may lead to wrong result in the analysis of the gauge 
theory. Here we define the ladder approximation in the SD equation (for the fermion 
propagator) as an approximation in which the full vertex function is replaced with 
the bare one. In our definition of ladder approximation, the gauge boson propagator 
is not necessarily restricted to the bare one. 

However, in our approach we look for a set of consistent solutions for A(p 2 ) and 
B{p 2 ) together with the vertex function T^(p,q) and the gauge boson propagator 
D/j.v(k) under the adopted approximation, so that the solutions are consistent with 
the WT identities. For the bare vertex approximation to be consistent with the WT 
identity ( 2.16j ), therefore, there should be no wavefunction renormalization for the 



fermion. In other words, the bare vertex approximation should self-consistently yield 
the result: 

A{p 2 ) = 1, (2.19) 



2 This naming is somewhat misleading, since the SD equation for the fermion propagator does 
not have an exact graphical representation corresponding to ladders. 
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as a solution of SD equation Eq. (|273|). In QED with covariant gauge ( |2.7| ), ( p,19| ) 



is realized only if we take both the Landau gauge and the quenched ladder approx- 
imation where the full photon propagator is also replaced with the bare one, i.e. 
I1(A;) = 0. [] Therefore, within the framework of SD equations in QED, the bare 
photon propagator in the Landau gauge together with the bare vertex is a set of 
consistent solutions of the SD equation (the fermion propagator can be non-trivial). 
[] This fact is well-known in four dimensions ]27| and holds in any dimension D > 2, 

Hi- 



sec 



If gauges other than the Landau gauge are adopted, A(p 2 ) = 1 does not hold even 
in the quenched ladder (bare vertex and bare gauge boson propagator) approximation 
and we have to solve the coupled equations for A and B. However, it is known that 
the solutions in QED obtained in such a scheme are severely gauge-dependent, see e.g. 
p9| . Indeed, the simplest quenched ladder approximation does not satisfy the WT 



identity except for the Landau gauge in the sense described above. In order to obtain 
gauge-parameter independent results in this case, we need to modify the full vertex 
so as to satisfy the WT identity. In the framework of the SD equation, however, it 
is rather difficult to obtain gauge-parameter independent results by modifying the 
vertex function. [] 

As explained above, in the quenched ladder approximation, only the Landau gauge 
t](k 2 ) = 1 can give a set of consistent solutions for the SD equation in this sense. How- 
ever, if the vacuum polarization is included in the photon propagator (i.e. unquenched 
case n(/c) ^ 0), the Landau gauge has no longer this property in the SD framework. 
Adopting the nonlocal gauge enables us to extend this scheme beyond the quenched 
ladder approximation. In a gauge theory, there is the freedom of choosing such a 
gauge. The existence of a gauge where A(p 2 ) = 1 in the SD framework was shown 



independently by Georgi, Simmons and Cohen [17] and by Kugo and Mit chard 



in four-dimensional gauge theory. The extension to arbitrary dimension D(> 2) is 



straightforward, as done in the Appendix of plfl . There the bare vertex was assumed 
from the beginning. In the next section we re-derive the nonlocal gauge in a more 
general setting. 

As long as the full vertex has the form ( p.!4| ), the nonlocal gauge plays the same 
role as the Landau gauge of the quenched ladder QED. In the next section, we show 
that the function r)(k 2 ) can be chosen so that the SD equation Eq. ( |2.3|) for the fermion 
propagator with the vertex function ( j2.14 ) leads to the solution Eq. (|2.19| ) for A(p 2 ), 



i.e. no wavefunction renormalization. Therefore, the nonlocal gauge gives the most 
economical choice, since we have only to solve the single equation for B. 



3 In the usual perturbation theory of QED, the Landau gauge is a special gauge in which the 
vertex correction vanishes. The above example shows that the similar situation occurs also in the 
framework of the SD equation. 

4 It is argued (30) that the quenched limit of gauge-invariant study of QED4 coincides with 
the Landau-gauge results obtained from the Schwinger-Dyson equation in the quenched ladder (bare 
vertex) approximation. This suggest that the Landau gauge is the best one in this approximation. 

5 In order for the full vertex function to be consistent with the Ward-Takahashi identity, the 
ansatz for the vertex function should include S, e.g. r^(p, q) = ^■[S^ 1 (p) — S^ 1 (q^\ + Tj t (p, q) . 
However, this requirement is not sufficient to determine the vertex uniquely, see ref. |32[ |33| , p4f and 
references therein. For an other proposal, see [p5[. 
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So far, we have not discussed the explicit dependence of G on A. Examples used 
so far for the ansatz belonging to the class ( [2.1 4|) are as follows. 



G(p 2 ,q 2 ,k 2 ) = A(p 2 ), A(q 2 ), 

= \W)+A{q% 



A{p 2 )9{p 2 -q 2 )+A{q 2 )0{q 2 -p 2 ), 

A{p 2 ) + Cl p 2 A'{p 2 ) + c 2 p A A"(p 2 ) + ...} 6{p 2 - q 2 ) + {p^ gX2.20) 



with some constants q(i = 1, 2, ...). If we take the nonlocal gauge so that A{p 2 ) = 1, 
the function G(p 2 , q 2 ,k 2 ) reduces to 1 uniformly in p, q, k (when B{p 2 ) is neglected), in 
order to be consistent with the WT identity. The following ansatzes are incompatible 
with the WT identity: 



G(p 2 ,q 2 ,k' 



A n (p 2 ), A n (q 2 ), 

^{A(p 2 ) + A(q 2 )T, 

A(p 2 )A(q 2 ), 
A(p 2 )A(q 2 )/A(k 2 ), 



(2.21) 



for any integer n > 2. 

In the massive fermion phase B(p) ^ 0, the small deviation of G(p 2 , q 2 , k 2 ) from 1 
is O(B) which gives rise to at least 0(B 2 ) terms in the integrand of the SD equation 
for B. In order to study the critical behavior in the neighborhood of the critical point 
through the solution of the SD integral equation for B(p 2 ), we can put G(p 2 , q 2 , k 2 ) 
I in the integrand of the integral equation after taking the nonlocal gauge, even 
if G(p 2 , q 2 , k 2 ) might deviate from 1 off the critical point. This is because the order 
0(B 2 ) terms are irrelevant to the bifurcation solution (from the trivial one B(p 2 ) = 0) 
which is sufficient to study the critical behavior. 

Therefore, in this approach we do not need the explicit form of the vertex in 
order to study the critical behavior of the model. This is a further advantage of this 



approach. For other types of the full vertex with different tensor structure |32|, |33|, |34 



such a convenient gauge is not known and we must solve the coupled equation for A 
and B as well as T^. In our approach all effects coming from the vertex correction are 
incorporated into the SD equation by modifying the gauge boson propagator through 
the nonlocal gauge. 

In QED3 it has been shown that the nonlocal gauge reproduces systematically 
the previous result obtained for the fermion self-energy with corrections up to the 
next-to-leading order in 1/N expansion 0, but in the approach of |J it was necessary 
to consider the vertex correction in writing down the SD equation. 
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3 Derivation of nonlocal gauge 

We decompose the SD equation into a pair of integral equations according to the 
following procedure: 



Atf) = 1 + ^ 



B(p 2 ) = 

where S denotes the self-energy part: 

d D q 



p 2 tr(l) 
tr[S(p 2 )] 



S (P 2 ) j^^lf,S{q)T u (q,p)D^{p~q). 



(3.1) 
(3.2) 

(3.3) 



We use the ansatz ( 2.14 ) for the vertex function. 

Then the SD equation (|3.1| ) for the fermion wave function renormalization A reads 



p 2 A[p 2 ) -p 2 

d D g A(q 2 )G(p 2 ,q 2 ,k 2 ) 
(2n) D q 2 A 2 (q 2 ) + B 2 (q T 



xk 2 D T (k 2 



(D-2] 



p-q fp-q p 2 q 2 - (p ■ q) 



k 2 



+ 



k 2 



- 2- 



k- 



V (k 2 



(3.4) 



On the other hand, the SD equation (|3.2|) for the fermion mass function B reads 

d D q B(q 2 )G(p 2 ,q 2 ,k 2 ^ 



B( P 2 



m + e 



[2n) D q 2 A 2 {q 2 ) + B 2 {q 2 
Separating the angle $ defined by 



-D T (k )[D - T](k )]. (3.5) 



k 2 := (q — p) 2 = x + y — 2^/xy cos??, x := p 2 , y := q 



(3.6) 



we find (for D > 2) 

xA(x) — x 

fA2 y (D-2)/2 A ^ 



C D e 2 I dy 
Jo 

xk 2 D T (k 2 ) 



yA 2 (y)+B 2 (y) 



I* d$sm D ~ 2 $G(p 2 ,q 2 ,k 2 ' 
Jo 



cos$y/xy{D — 2 + rj{k 2 )] xy — (y/xycosfi) 2 



k 2 



k 1 



where 



C D :~- 



2r» 7r (r»+i)/2r(^l > 



V (k 2 ) , (3.7) 



(31 



We follow the same procedure as that given in [|18j and Appendix of ref. [21]. We 
perform the angular integration by parts according to 



dfisin cos i^y/xyf (a 



xy 



D - 1 



sm D ~ l $f(z) 



7T 2xy 
o ~ D - 1 Jo 







d$sm u $—f(z) 
oz 



(3.9) 
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where z := k 2 = (q — p) 2 and the differential with respect to z is done with x and y 
being fixed. Thus we find 

xA(x) — x 

d§ sin D d 



C D e 2 /-A 2 y^- 2 )l 2 A{y){2xy) 
dy- 



D-Uo ~* yA\y) + B\y) 
x ^ y, z) [D - 2 + 7j(*)]£> T (*)} + (£) — ! 



G(x,y,z)D T (z)r)(z) 



(3.10) 



This is further rewritten as 



xA(x) 



x 



C D e 2 /-a 2 dy y( D - 2 y 2 A(y)(2xy) 



D — 1 Jo 
1 



y D-l 



G(x,y,z) j 

dG(x, y, z 
dz 



yA\y) + B\y) Jo 

)(- 

Mas 



citf sin 13 



■[D-2 + r/(^)]z D - 1 J D T (z 



From this equation, it turns out that the requirement A(p 2 ) 
tive of 5, if 



dG(x,y,z) 
dz 







(3.11) 

1 is achieved irrespec- 
(3.12) 



(3.13) 



and 7](k 2 ) satisfies the following differential equation: 

-^-[z^DrizHz)] = —(D - 2)z D ~ 1 ^. 

Thus, once the function Dj<{k 2 ) is given, we can find the nonlocal gauge 7](k 2 ) by 
solving Eq. (|3.13|) , so that A(k 2 ) = 1 follows under the ansatz for the vertex function 
( |2.14| ), provided that G does not depend explicitly on k 2 . Then we have only to solve 
Eq. (|3.5| ) for the fermion mass function B(p 2 ). 

d D q B(q 2 ) 



B(p 2 ) 



m + e 



or 



f 

B(x) = m + Cdg 2 / dy 
Jo 

where the kernel is given by 



(2n) D q 2 + B 2 (q 2 ) 

A 2 , y^l 2 B(y) 



y + B 2 (y) 



G(p 2 ,q 2 ,k 2 )D T (k 2 )[D- V (k% (3.14) 
G(x,y,z)K(x,y), (3.15) 



K(x,y):= / d$ sin u - 2 $D T (k 2 )[D -7i(k% 



(3.16) 



The quenched ladder QED in the covariant gauge is recovered by taking U(k) = 0, 
i.e. Dt(z) — 1/z. In this case the nonlocal gauge reduces to the local gauge: rj(k 2 ) = 
i] = 1 — £. This reproduces the well-known result: A(p 2 ) = 1 in the gauge rj = 1, i.e. 
Landau gauge £ = for D > 2. Especially, for D = 2, A(x) = 1 is satisfied in the 
gauge 77 = 0, i.e. Feynman gauge £ = 1, as can be seen from Eq. (|3.4j ). However, in 
D = 2, the SD equation can be exactly solved in an arbitrary gauge, see [f2~5|. 
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3.1 Nonlocal gauge without IR cutoff 

The differential equation ( |3.13| ) is a first order differential equation and is simply 
solved by choosing a boundary condition. Integrating both sides of Eq. (|3.13| ) from 
to k 2 , we obtain 

where the prime denotes the differential with respect to z. Here we have assumed the 
boundary condition: 

[z D - 1 D T (z) V (z)]\ z=0 = 0, (3.18) 
so as to eliminate the l/z D ~ l singularity in T](z). Alternatively, we can write 

D - 1 



V (z) = (D-2) 
where we have assumed that 



z D - l D T (z) Jo 



f dtD T (t)t D - 2 -l 
Jo 



(3.19) 



[z D - 1 D T (z)}\ z=0 = 0. (3.20) 
This should be checked after having obtained the function r](z). 

3.2 Nonlocal gauge with IR cutoff 

For later convenience, we introduce an IR cutoff e in the nonlocal gauge by integrating 
fl3jj) from e 2 to k 2 : 

v(k 2 ) - ^^K U 2 ) D ~ 2 f dzP> (z)z D -i (3 21) 

V{k) ~ (k 2 ) £> - 1 D T (k 2 ) (k 2 ) D -iD T (k 2 ) U dZ r{Z)Z ' {S } 

where t](e 2 ) is undetermined. Note that Eq. (|3.13|) is rewritten as 

( , z^DAzWjz) + (D — 2)z D ^D' T {z) 

m ~ [z D -w T ( z )y • {6 - 22} 

Therefore, if we impose the flatness condition i]'(e 2 ) = on t](k 2 ) at k 2 = e 2 as a 
boundary condition, []the value r/(e 2 ) is determined as 

, 2 , (D-2)(e 2 ) D - 1 D'(e 2 ) 

Hence we arrive at the expression of the nonlocal gauge: 

(D - 2)e 2 ^D T (e 2 ) (e 2 ) " 1 D' T (e 2 ) 



7]{k 2 



(k 2 ) D -w T (k 2 ) [z D -w T (z)}'\ z=£ i 

/ dzD' T {z)z D -\ (3.24) 



(k^-Wrik 2 ) Je* 



This means the flatness of the effective coupling, see section 6. 
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or 

e 2(£-l)£) r ( e 2) - i,-2\L>-\. TV !/2, 



v (k 2 ) = 1 - ^) = -(D-2) 



{k 2 ) D -W T (k 2 ) 



[z D -W T (z)}'\ 



D~2)+ %J^-r^ 1^ dzD T {z)z D ~ 2 . (3.25) 



(k 2 ) D -W T (k 2 ) 

Note that £(fc 2 ) := 1 - r](k 2 ) is in general different from £(fc 2 ). 

4 Gauge choice and LK transformation 

It is worth remarking that the SD equation and the WT identity are form-invariant 
under the the Landau-Khalatnikov (LK) transformation fi~6|| : 



D'^(x) = D, u (x)+d,d u f(x), 
S\x,y) = e e2[f{o) - f{x - y)] S(x,y), 
V v (x,y,z) = e e2 ^-H x -^V„(x,y,z) 

+S(x, y)e e2 ^-^-^d z u [f(x -z)- f(z - y)}, (4.26) 



where 



D^{x,y) = (0\T[A,(x)A u (y)}\0), 

S{x,y) = (0|r[*(a;)%)]|0), 
V„(x,y,z) = (Q\T[*(x)*(y)A v (z)]\0). (4.27) 

This can be easily shown in coordinate space where the SD equation has the following 
form: 



and 



(id - m)S{x, y) = 5 D (x - y) + ie 2 ^(${x)^{y)AJx)}, (4.28) 



D~Xx,z) = D^- 1 (x ) z)-U^(x ) z) 



n^(x, z) = (g^ud 2 - d^d v )Yi(x - z) 



ie 2 



J d D z 1 d D z 2 ti[^S(x, z 1 )T v (z 1 , z 2 ; z)S(z 2 , x)}, (4.29) 
where 

(^(x)^{y)A fl (z)) = f d D x'd D y'd D z'S(x,x')T u (x',y'-z')S(y',y)D, u (z',z). (4.30) 



Therefore, if we know a consistent set of solutions (for the full gauge boson propagator, 
the full fermion propagator and the full vertex function) of SD equation in a single 
gauge, the solutions in other gauges are obtained through the LK transformation. 

So far, the LK transformation has been used to transform the Landau gauge 
result £ = into other gauges £ ^ 0. Note that it is possible to perform the inverse 
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LK transformation (from non-Landau gauge to Landau gauge). It turns out that 
the inverse LK transformation is obtained from ( [4.26 ) by replacing f(x) with —f(x). 



Furthermore, the LK transformation and its inverse allows us to deal with the nonlocal 
gauge fixing, since f(x) is an arbitrary function. Therefore, the Landau gauge result 
is recovered from the nonlocal gauge by choosing 

where £(fc 2 ) is related to 7](k 2 ) by (|2T3l) . The inverse LK transformation enables us 
to compare the result of this approach with the conventional one. 

Using the inverse LK transformation, we can obtain the fermion propagator Sl(x) 
in the (usual) Landau gauge £ = from the fermion propagator S n i g (x) in the nonlocal 
gauge 

S L (x) = e A ^S nlg (x), (4.32) 

where 

= -WW. /<*> - w=Bm- (4 ' 33 > 

In general, the fermion propagator in configuration space is written in the form: 

S(x) = i 7 %P(x) + Q(x), (4.34) 
in accord with S(p) = [A{p 2 )^p^ + B(p 2 )]" 1 . Taking into account the relation 



we obtain 



B(p 



2\ 



A 2 (p 2 )p 2 + B 2 (p 2 ) 



J d D xe ip - x Q{x). (4.36) 



Therefore, A(p 2 ) and B{p 2 ) in the Landau gauge are obtained by substituting P(x) = 
e M x ) p rdg (x) and Q(x) = e A ^Q n i g (x) into ( |4.36| ). In order to accomplish this, we must 



obtain the expression of the fermion propagator S n i g {x) in configuration space. 

In the massless fermion phase B(x) = 0, the fermion propagator S n i g (p) in the non- 
local gauge is nothing but the free massless propagator S n i g (p) = So(p) = l/(7 M p M ). 
In configuration space the free massless propagator is given by 

§o{x) = / Wf^'^ = ^mW, Port = Ji/^ /2 - (4-37) 
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Therefore, the full fermion propagator in Landau gauge is given by 

S L (p) = ^bt, A ~\P 2 ) = * J d D xe^ x e A ^\p ■ x)P (x). (4.38) 

Note that the vacuum polarization IT(fc) determines £(k) and then A(x). Therefore, 
the specification of the vacuum polarization is crucial. In the quenched limit H(k) = 0, 
= and hence Sl = So- This is a consistency check. 
In the massless fermion phase, the vacuum polarization can be calculated easier 
than the massive case. The explicit expression depends crucially on the spacetime 
dimension in question. The nonlocal gauge function is expected to have a finite 
range, namely, \£(h)\ < c uniformly in k, except for the neighborhood of k = 0. So 
n(/c) is expected not to change the UV behavior of f(k) qualitatively, see ( fp 



However, in lower dimensions 2 < D < 4, U.(k) can influence the IR behavior of f(k) 
considerably. This is discussed in section 5.3 in more detail. 

5 Running coupling constant and nonlocal gauge 
5.1 Wavefunction renormalization and running coupling 



Suppose that the SD equation (ft.4|) for A has been solved in the Landau gauge 
rj(k 2 ) = 1 under the ansatz ( |2.14| ) for the vertex function (when B 2 is neglected). 
Then the SD equation (|3.5| ) for B can be written as 



B(p 2 ) = m + e 2 J 



d D q B(q 2 ) 
(2tt) d q 2 



G(p 2 ,q 2 ,k 2 ^ 



A 2 (q 2 



(D-l)D T (k 2 ), (5.1) 



where B 2 (q 2 ) in the denominator is neglected according to the bifurcation method 

H p. 

In the paper the following ansatz for the vertex is adopted in order to look 



for the solution in QED 3 : 

G(p 2 ,q 2 ,k 2 ) = A(q 2 r, (5.2) 

where n is an integer (n = 1,2, ...). Under this ansatz, the approximate solution of 
the SD equation for A was obtained in the usual Landau gauge: 

A(p 2 ) = (l + ~^Kt) , (5.3) 

where 



7T z iVf 



and 



t:=ln^, a:=^£. (5.5) 
a 8 
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Based on this solution, it was pointed out that one can define a running coupling 
constant K(t) whose actual running is given by 

K(t) := A:f(p*)K = 1 + ^ Kf (5.6) 

due to wavefunction renormalization. Indeed, in the absence of wavefunction renor- 
malization Aijp 2 ) = 1 and the coupling K(t) does not run, i.e. K(t) = K. 

It is claimed that the running coupling constant of QED 3 obtained in such a 
way corresponds to the asymptotically free case for n < 2 (here n — 1 is likely to 
be the physical one due to the WT identity). If we accept Eq. (|5.6|) at face value, 
the running coupling K(t) diverges in the IR region for n < 2. In this case, K(t) 
becomes strong enough to be able to cause chiral symmetry breaking and to make 
the bound state It should be noted that the asymptotic freedom defined here is 
valid only if t := lne/a < := — 3/[(2 — n)K). Therefore this asymptotic freedom 
would disappear for sufficiently large IR cutoff e or small UV cutoff a. This gives a 



possible (physical) explanation for the controversy among the results |14|, [8] on the 



phase structure of QED 3 . For more details, see ref. jll 

The above observation is based on the approximate solution ( |5.3| ) for A in the 
Landau gauge £ = under the vertex ansatz (|5.2j ). Quite recently, the rate of run- 



ning of the coupling constant K{t) has been studied in more detail by Aitchison and 



Mavromatos p3| and a subsequent paper |2J| where the validity of all the approx- 



imations made in solving the SD equation for A in the paper [TJJ were reexamined 
thoroughly under the same vertex ansatz. 

5.2 Running coupling through nonlocal gauge 

Now we study the running coupling based on the nonlocal gauge. In the nonlocal 
gauge, the function B obeys the SD equation ( p,14j ). If we choose the nonlocal 



gauge £(fc) such that A{p 2 ) = 1, the function G(p 2 ,q 2 ,k 2 ) should be replaced with 
1, i.e. G(p 2 ,q 2 ,k 2 ) = 1 from the consistency with the WT identity. In this case, the 
wavefunction renormalization disappears and the SD equation for B reduces to 

B(P 2 ) = m + e 2 / gL _ B(g V> ~ ^)}D T (k\ (5.7) 



This should be compared with the SD equation ( |5.1|) . Then the running coupling 
constant is given by 

miK = = £^i±m. (B . 8) 

Here note that the argument of the running coupling constant is the gauge-boson 
momentum, which is required for the nonlocal gauge to be consistent with the axial 
WT identity as well as the vector WT identity [fT8 1 . Therefore, the problem of finding 



the running coupling constant reduces to finding the nonlocal gauge. 

Therefore we can see the following advantage or benefits of this approach: 
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We do not have to assume a particular ansatz for the vertex. We can study a 
class of vertices of the form: r M (p, q) = 7 At G(p 2 , q 2 , k 2 ), without assuming any- 
particular form of G(p 2 , q 2 , k 2 ). The consistency with the WT identity requires 
that G(p 2 ,q 2 ,k 2 ) ^ 1 as A(p 2 ) — > 1. This class of vertex ansatz includes the 
previous one (|5.2| ) with n = 1. 



2) Under this ansatz for the vertex, we need not to solve the SD equation for 
A(p 2 ). This releases us from worrying about the validity of a number of ap- 
proximations which are required to solve the SD equation for A. Instead, we 
obtain the differential equation for the nonlocal gauge. This is solved by simple 
quadrature. In particular, we do not have to perform any angular integration. 



So the Higashijima-like approximation [38] is unnecessary to separate the kernel 



of the integral equation (at least in the normal phase). 

3) If necessary, we can recover the solution in the usual gauge (e.g. Landau gauge), 
by making use of the (inverse) LK transformation as discussed in the previous 
section. This enables us to compare the result in this approach with the con- 
ventional result. 

4) We can study the effect of a cutoff in more detail. 

The final point needs more explanations. In this approach it is not necessary to 
introduce the UV cutoff a in QED 3 . In the analysis of QED 3 the existence of 



the IR cutoff e was essential as well as the UV cutoff a. Whether or not there exists 
a "finite" critical number of flavors Nj above which (Nf > Nj) spontaneous chiral 
symmetry breaking disappear depends crucially on the existence of the infrared cutoff 
and the ratio e/a. It is expected that the asymptotic freedom claimed above would 
disappear for sufficiently large IR cutoff e or small UV cutoff a and hence this leads to 
a finite critical number of flavors iVj < oo for sufficiently large e/a. This observation 
is based on a naive analogy with QCD in 3+1 dimensions (QCD4). In QCD4, the 
theory has only one phase, the chiral-symmetry breaking and confining phase, due 
to asymptotic freedom where the running of the gauge coupling constant is essential 
On the other hand, QED 4 is not asymptotically free and has a critical gauge 



coupling constant e c above which (e > e c ) a fermion mass is dynamically generated 
and chiral symmetry is spontaneously broken. Therefore, the IR behavior of the 
running coupling constant in QED3 is very important to resolve the phase structure 
of QED 3 . 

On the the hand, QED 3 has already been studied in the nonlocal gauge [21], ^] 
where we have found a phase transition at a finite critical number of flavors Nj < 00 
which separates the chiral-symmetric phase from the spontaneous-chiral-symmetry- 
breaking phase. In this analysis we did not introduce the infrared cutoff from the 
beginning. Apparently, the results of [JTTJ] and [^l[ contradict with each other, since 



the very small IR cutoff (compared with UV cutoff a) should lead to the asymptotic 
freedom and no phase transition, i.e. non-existence of a finite critical number of 
flavors according to [fLl |. In this paper we resolve this apparent contradiction. For 



this, we introduce the infrared cutoff in the setting up of nonlocal gauge. The result 
is given in the next sections. 



16 



5.3 Inverse LK transformation and wavefunction renormal- 
ization 

Thanks to the inverse LK transformation, we can obtain the fermion propagator Sl(x) 
in the (usual) Landau gauge £ = from the fermion propagator S n i g (x) in the nonlocal 
gauge. In D = 3 dimensions, the free massless propagator in the configuration space 
is given by 

S (x) = (5.9) 

For D = 3, the massless fermion generates the vacuum polarization (at one- loop): 

Il(Jfe) = -ak. (5.10) 

By substituting the vacuum polarization function ( |5.1(J| ) and S n i g (x) = S (x) 
given by ( |5.9| ) into (|4.32| ) with the nonlocal gauge function £(/c 2 ) obtained explicitly 
in the next section, we can get the fermion propagator Sl{x) in the Landau gauge in 
configuration space. 

If £ was a constant, we would have for D = 3 

A(x) = e 2 / ^(e*- - 1)4 = (5-11) 



(2tt) 3V 



57T 



According to ([4.38 ) , it is not difficult to show that the wavefunction renormalization 
function is obtained as 

A-\ V 2 ) = 1 - ^l a rctan(4?). (5.12) 



In the IR region, we find that 

A-\p 2 )= l -{^? + 0^). (5.13) 
3 



This is totally different from what we expect based on the nonlocal gauge. 

In the nonlocal gauge, the integrand f(k) of A(x) in the UV region k — > oo 
exhibits the behavior 

/(*)~^T->0 (*-oo), (5.14) 

so the effect of LT is neglected in this region. However, in the IR region k —> 0, the 
presence of H(k) totally changes the situation: 

(k^O). (5.15) 

ak 6 

In either case, the situation does not resemble the constant £ case. If the integra- 
tion is performed numerically in (|4.38|) , we will be able to obtain the wavefunction 
renormalization function A(p 2 ) in the Landau gauge (the corresponding vertex func- 
tion can be also obtained from the LK transformation). The numerical result will be 
given elsewhere. 
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6 Running coupling constant of QED3 

In D = 3 dimensions, the running coupling (|5.8|) is obtained by shifting and scaling 
the nonlocal gauge according to 

K(t)/K = l + $yp-. (6.1) 

The exact non-trivial wavefunction renormalization in the Landau gauge is obtained 
from ( (4.38|) . Under a specific ansatz (|5.2| ), it is obtained from (|5.6|) and ( [5.8D as 



^ 2 )=(l + ^) 1/(n \ (6.2) 
and in particular for n = 1 

A(P 2 ) = (l + ^) 1 • (6-3) 
Therefore, we study the behavior of the nonlocal gauge function £ (k 2 ) in what follows. 

6.1 Running coupling constant (I) 

From Eq. ( |3.25| ), the nonlocal gauge of QED3 with an IR cutoff e is given by 



(k 2 ) 2 D T (k 2 
e(k 2, ._ ^D T (e 2 ) \ (e 2 ) 2 D' T (e 2 ^ 



(k 2 ) 2 D T (k 2 ) 



[z 2 D T {z)]'\ z -_ 



- 1 



(6.4) 



Here we decomposed the nonlocal gauge into two pieces: the first piece reduces in 
the limit e — > to the nonlocal gauge without an IR cutoff, while the second piece £ b 
comes from the flatness condition £'(e 2 ) = and vanishes in the limit e — > 0. First of 
all, we consider the following integral 



Ji(k 2 ;e) := / dz 



* 2 z 



2; + a-y/i 

k 2 -2aVk 2 -e 2 + 2aV?+2a 2 \n - . (6.5) 

V e 2 + a 



This has the expansion: 



JJk 2 ; e) = (-e 2 + 2 a v 7 ? + 2 a 2 log(a) - 2 a 2 logfa + v^ 2 )) + — + 0(A; 4 ). (6.6) 

3 a 

Hence the first piece £ a of the nonlocal gauge (|6.4j) given by 



C(fc 2 )=2-2 fc2 + ^ J 1 (A: 2 ;e 2 ) (6.7) 
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has the following expansion around k = 0: 

2 a (e 2 - 2aV^+2a 2 log(l + \fe*/a) 



k 3 

2 (e 2 - 2 a + 2 a 2 log(l + yf^/a) 



k 2 

2 k k 2 k' , „ a „ , 

+ o -T^ + ^-T^ + O^ 4 ). (6i 



3 3 a ' 5 a 2 15 a 3 
Note that the introduction of the IR cutoff e generates IR singular terms like 2ack~ 3 , 
2ck~ 2 where the coefficient c is always positive. Such a singular behavior in the IR 
region disappears if we put e = from the beginning, and £ a reduces to 



2 k k 2 2 k 3 

3 3 a 5 a 2 15 a : 
This is nothing but the result obtained in P] (by setting the Chern-Simons coefficient 



9 equal to zero: 8 = in eq. (29) of |2T|). On the other hand, the second piece is 
always negative and singular at k = 0: 

4 (k 2 + a VW) e 3 ^ / a 1 \ 

i b e {k 2 ) = K Lr- = -7 —=r ^ + ^ ) . (6.10) 

fc 4 (3a + 2v^) (3a + 2v / ?) V^ 3 k ^ 

By adding (|6.^ ) and ( 6.10 ), we find that the singular part with negative power of k 
in the nonlocal gauge is negative. 

In the region k/a 3> 1, £ a is dominant, because £ b decreases more rapidly than £ a 
which behaves as 

^(A; 2 ) = 2a | 4a 2 + 26 2 - 4av^ + 4a 2 log((a + y^)/A:) | 1 3 n 
k k 2 k 

It turns out that the alone is monotonically decreasing in k with a maximum value 
2 at k = e (and diverges monotonically as k — > 0) which rather enhances the effective 
coupling compared with £o(0) = 2/3. This tendency agrees with the observation 
made in the previous paper |TT[ where the running of the coupling is terminated at 
k = e and the flatness of the coupling K(k 2 ) = K(e 2 ) for k < e is assumed a priori, 
which is borrowed from the QCD 4 analysis. If we consider only the £ a piece, there is 
a discontinuity in the derivative £,f{k) at k — e, since + < and — = 0. 
However, the inclusion of 0) which is necessary to satisfy ^(e) = (continuity 
of £' e (k) at k = e) considerably changes the situation. The inclusion of such a term 
causes a slowing down of the rate of decrease of the effective coupling constant K{t). 



This tendency agrees with the recent analysis of [24]. The the nonlocal gauge function 



and the running coupling constant are monotonically decreasing in k for k > e and 
have upper bonds: 

Ue 2 ) = ^(e 2 ) + £ e V) = 2 - = — %-r- < I (6.12) 

3a + 2e 3 + 2e/a 3 

For k < e, the nonlocal gauge £ e decreases as k decreases. In Figure 1, the nonlocal 
gauge is plotted for e/a = 0.1. The running coupling constant is obtained from ( |6.1|) , 
and the wavefunction renormalization in the Landau gauge is obtained from (|6.3j). 
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6.2 Running coupling constant (II) 

In this section, we consider another way of introducing the IR cutoff. We introduce 
the IR cutoff 5 in the gauge-boson propagator: 

- wrhrp- (6 - 13) 

The cutoff 5 plays the same role as the gauge-boson mass and seems to be more 
natural than the previous one. This choice of Dt is equivalent to the gauged Thirring 
model in the nonlocal gauge [H5 where 5 2 = e 2 G^ 1 for the Thirring coupling Gt- 
In this case, the nonlocal gauge reads 

is(k 2 ) = 2 - 2 k2+ ^ + P J 2 {k\ 5), (6.14) 

with 

rk 2 



J 2 (k 2 ;5) := / dz 
Jo 



Z + CXy/z + 5 2 



k 2r 3 
dr- 



o r 2 + ar + 5 2 
(k - 2a)k + (a 2 - 5 2 ) \n\(k 2 + ak + 5 2 )/5 2 \ 

+a(36 2 - a 2 ) \I B {k) - I B {0)), (6.15) 



where the indefinite integral is(/c) is defined by 

I B (k) := f k dr 1 (6.16) 
J r l + ar + o l 

The indefinite integral Isik) can be calculated: For a 2 > 4<5 2 , 

-i , (2t + a Wo 2 - tfy 



and for a 2 < AS 2 , 



Ib ^ = IaI 2 arctan = 2 • ( 6 - 18 ) 

V4cr — a 2 v4o — a 2 



When a 2 > 45 2 , a careful analysis shows that all the terms in J 2 up to 0(k 3 ) 
cancel, namely, 

■M* 2 ; i) = ^ - + + O(F). (6.19) 

This implies the following expansion of the nonlocal gauge £,&(k 2 ) around k = 0: 

&(*") = - 2a(5 ^- 6 1W2) ^ + O(^). (6.20) 
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This result shows that there is no singularity in the nonlocal gauge ( |6.14Q at k = 0, in 
sharp contrast with the nonlocal gauge £ e (k 2 ). In the region k/a <C 1, the nonlocal 
gauge behaves as follows: at 5 /a = 0.1 



l s (k 2 ) = 1 - 20(k/a) + 1300(fc/a) 2 - 92000(£;/a) 3 + 0((k/a) 



(6.21) 



and at 5 /a = 0.01 



= l - 2000(£;/a) + 1.333 x 10 7 (k/a) 2 - 9.52057 x 10 w (k/a 
+0((k/af). 



(6.22) 



It should be remarked that, in the limit (5 — > 0), the nonlocal gauge (|6.20| ) does 
not reduce to the expected form ( |6.9D . There is a discontinuity between the nonlocal 
gauge ^(fc 2 ) with an IR cutoff 5 and the nonlocal gauge £o(^ 2 ) obtained without an 
IR cutoff from the beginning. Such a discontinuity does not exist for the nonlocal 
gauge £,e(k 2 ) with IR cutoff e. Such a situation can be seen in the finite temperature 
case: there may occur some discontinuity between the zero-temperature limit of finite- 
temperature calculation and the corresponding result evaluated at zero-temperature. 
In this context 5 can be interpreted as a plasmon mass. This discontinuity is very 

121 



similar to the result found by Aitchison et al 
decreases more slowly than claimed in the previous paper 
is plotted in Figure 2. 



The non-local gauge £,s(k 2 ) 
The nonlocal gauge 
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6.3 RG-like point of view 

By using the idea of RG, we show that the IR and UV behaviors of the nonlocal gauge 
£,s{k 2 ) and the effective running coupling constant can be easily analyzed without 
performing any integration. For this, we calculate the derivative of the nonlocal 
gauge function and express the result in terms of the nonlocal gauge function itself: 

,,/;\ d ~ 9n 2k 2 + 3ak + 4<5 2 ~, 9n 

From this equation, we can observe the following: 

1) without IR cutoff {5 = 0), in the IR limit k -> 

P(£) = 4 + 3(i(k 2 ) - 2) = 3i(k 2 ) - 2. (6.24) 

This implies l{k 2 ) 2/3 = | (0) in the IR limit k ->• 0. 

2) with IR cutoff (5^0), in the IR limit k -»• 

(3(0 = 4 + 4(£> 2 ) - 2) = 4|(/c 2 ) - 4. (6.25) 

This implies £(& 2 ) — > 1 in the IR limit k — > 0. 

These result seems to show the existence of non-trivial IR fixed point for the 
running coupling constant K := K(t) defined by K(t)/K = l+l(k 2 )/2 at K = A/3K 
and K = 3/2K ~ 0(1/ N) corresponding to 1) and 2) respectively. 
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In particular, 

3) in the UV limit k — > oo (irrespective of the value of the IR cutoff), 

13$ = 4 + 2(£(k 2 ) - 2) = 2£> 2 ). (6.26) 

This implies £,{k 2 ) — > in the UV limit — > oo. 
Especially, 

4) in the quenched limit (Nf — > and 5 = 0), 

/?(£) = 4 + 2(£(* 2 ) - 2) = 2ft* 2 ). (6.27) 

for any k. This shows that the constant solution is possible: £(k 2 ) = 0. 

It is straightforward to generalize the above argument to arbitrary dimension. 



7 Conclusion and discussion 

In this paper we have discussed an alternative approach for obtaining the effective or 
running coupling constant and the RG property of gauge theory in the SD framework, 
which is appropriate for non-perturbative study of gauge theory. This approach can 
be applied only to a gauge theory, as we use the the gauge invariance as an essential 
ingredient in the approach. We argued that this approach is superior in several 
respects to the previous conventional approach in which the SD integral equation 
has been solved to obtain the wavefunction renormalization function A(p 2 ) under a 
specific vertex ansatz to define the running coupling. In this approach we do not solve 
the SD equation for A, instead we obtain the differential equation which the nonlocal 
gauge function must satisfy in order that A(p 2 ) = 1. 

The validity of this approach has been exemplified in the study of the IR behav- 
ior of the running coupling constant and non-trivial IR fixed point in QED3. This 
approach confirms a recent result: the slowing down of the rate of decrease of the 
running coupling constant and the existence of non-trivial IR fixed point, as claimed 
in [^3], [24]], but here without relying on a specific ansatz of vertex and on a number of 



approximations adopted to solve the SD equation in the conventional approach |[l 1 



In this paper we have not studied the dynamical mass generation and spontaneous 



chiral symmetry breaking |TJ] in QED3 by solving the SD equation for B in the 
nonlocal gauge. The nonlocal gauge £,&{k 2 ) does not qualitatively change the previous 
result |2T| in the nonlocal gauge £o(^ 2 ) without IR cutoff. For example, there should 



exist a certain finite critical number of flavors Nj above which dynamical mass is not 
generated and chiral symmetry is restored. 

The relationship between our approach and the conventional one can be seen by 
using the LK transformation. Under an LK transformation, the SD equation and 
the WT identity are form-invariant. By making use of an LK transformation, the 
non-trivial wavefunction renormalization A(p 2 ) in the Landau gauge can be obtained 
from the nonlocal gauge, as shown in section 4.3. Thus the conventional picture can 
be recovered by the (inverse) LK transformation from our approach, if desired. It 
would be interesting to apply this procedure to the case of QED3 for which Maris 
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I5| has recently given a rather full discussion of the SD equation in the conventional 



framework. 

The nonlocal gauge has been applied to gauge theories in four dimensions [|T7|, [TEfl , 

T9|j under a bare vertex approx- 



QED 3 |g, [23], |3 and the gauged Thirring model 
imation. Obviously, this approach is not so systematic as the perturbative method. 
The nonlocal gauge must be obtained case by case and there is no guarantee that 
such a gauge does exist beyond this order of truncation of SD equation. Neverthe- 
less it is sufficiently interesting to warrant the extension of this approach to finite 
temperature gauge field theory f3"9fl . In such a case, the nonlocal gauge may only be 
obtained approximately, as in the case for QED3 with a Chern-Simons term [El], |2"2"||. 
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Figure Captions 



Fig.l: Plot of nonlocal gauge given by ( |6.4| ) as a function of k/a. Three graphs corre- 
spond to i a {k) (above), ^{k) (middle) and £ e (/c) (below). Here we have chosen 
e/a = 0.1 

Fig.2: Plot of nonlocal gauge given by ( |6.14| ) as a function of k/a. Two graph corre- 
spond to e/a = 0.1 (above) and e/a = 0.01 (below). 
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